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Abstract 

We investigate analytically the problem of enumeration of nonequivalent prim- 
itive words in the braid groups Bn for n ^ 1 by analysing the random word 
statistics and target space on the basis of the locally free group approximation. 
We develop a "symbolic dynamics" method for exact word enumeration in locally 
free groups and bring arguments in support of the conjecture that the number of 
very long primitive words in the braid group is not sensitive to the precise local 
commutation relations. We consider the connection of these problems with the 
conventional random operator theory, localization phenomena, and statistics of 
systems with quenched disorder. Also we discuss the relation of the particular 
problems of random operator theory to the theory of modular functions. 
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1 Introduction: Problems and Motivations 



Recent years have been marked by the emergence of more and more problems related 
to the consideration of physical processes on noncommutative groups. In trying to 
classify such problems, we distinguish between the following categories in which the 
noncommutative origin of phenomena appear with perfect clarity: 

1. Problems connected with the spectral properties of the Harper-Hofstadter equa- 
tion dealing with the electron dynamics on the lattice in a constant magnetic field. We 
mean primarily the consideration of groups of magnetic translations 0, § and properties 
of quantum planes 0. 
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2. Problems of classical and quantum chaos on hyperbolic manifolds: spectral proper- 
ties of dynamical systems and derivation of trace formulae ^, ^ as well as construction 
of probability measures for random walks on modular groups [Q. 

3. Problems giving rise to application of quantum group theory in physics: deforma- 
tions of classical abelian objects such as harmonic oscillators |^ and standard random 
walks ITD . 



4. Problems of knot theory and statistical topology: construction of nonabelian 
topological invariants |jll], |12|, consideration of probabilistic behavior of the words on 
the simplest noncommutative groups related to topology (such as braid groups) |ll3| , 



statistical properties of " anyonic" systems 



5. Classical problems of random matrix and random operator^ theory and localiza- 
tion phenomena: determination of Lyapunov exponents for products of random non- 
commutative matrices [1^, [l^, study of the spectral properties and calculation of 
the density of states of large random matrices []TB|, . 

Certainly, such a division of problems into these categories is very speculative and 
reflects to a marked degree the authors' personal point of view. However, we believe 
that the enumerated items reflect, at least partially, the currently growing interest in 
theoretical physics of the ideas of noncommutative analysis. 

Let us stress that we do not touch upon the pure mathematical aspects of noncom- 
mutative analysis in this paper and the problems discussed in the present work mainly 
concern the points 4 and 5 of the list above. 

It is widely considered that a new fresh stream in topology was brought about by 
the recognition of the fact that there exists a deep relation between the Temperley-Lieb 
algebra and the Hecke algebra representation of the braid group. This fact resulted in 
the remarkable geometrical analogy between the Yang-Baxter equations, appearing as 
a necessary condition of the transfer matrix commutativity in the theory of integrable 
systems on the one hand, and the Reidemeister moves, used in the knot invariant con- 
struction on the other hand. We can mention several nice reviews [IT, 2T| and books 

^Following L.A. Pastur, we will distinguish random matrices and matrix representations of random 
operators. To the random operators we attribute the n x n tables having of order n random entries; 
if the number of random entries grows faster than n when n — > c», we call such a table as a random 
matrix. 
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22| , [T^ regarding the construction of new knot and link invariants in terms of inte- 
grable 2D statistical models, as well as their relation with different matrix and tensor 
representations of some noncommutative groups. 

Besides the traditional topological issues such as construction of topological invari- 
ants, investigation of homotopic classes and fibre bundles, we consider a set of adjacent 
but much less studied problems lying at the border between statistical physics, probabil- 
ity theory and topology. First of all, we should mention problems related to the so-called 
"knot entropy" calculation. Most generally this set of problems can be formulated as 
follows. Take the lattice embedded in the space R^. Let f^Ar be the ensemble of all 
possible closed non-self- intersecting A^-step loops with one common fixed point onZ^; 
by uj we denote the particular trajectory configuration. The question is: what is the 
probability Vn of the trajectory uj eVL]^ belonging to some specific homotopy class. In 



13| it has been shown that many non-trivial properties of the limit behavior of knot 



statistics can be explained in the context of the limit behavior of random walks over the 
elements of some nonabelian (hyperbolic) group related to the braid representation of 
knots. 

In the context of "topologically-probabilistic" consideration, the problems dealing 
with the limit distributions of noncommutative random walks were not discussed prac- 
tically except very few specific cases |]16|, 0, |2^, In particular, in these works it has 
been shown that statistics of random walks with a fixed topological state with respect to 
the regular arrays of obstacles on the plane can be mapped to random walks on the free 
group r2 which has the topology of a simply connected tree. The analytic construction 
of nonabelian topological invariants for trajectories on a double punctured plane as well 
as the statistics of the simplest nontrivial random braid i?3, were discussed briefly in 



25, ra. 



A preliminary analytical and numerical study of the statistics of random walks 
(Markov chains) on braid and so-called "locally free" groups | (see definition below) 



was recently undertaken in works In the case of the braid group, the rather 

complicated group structure prevents us from applying the simple geometrical pictures 
of the free group r2 (see PPI). Nevertheless the problem of the limit distribution for 



random walks on can be reduced to the problem of a random walk on some graph 



27, 28]. In case of the group i?3 we were able to construct this graph explicitly, whereas 



^This notation has been introduced by A.M. Vershik in 
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for the group Bn {n > 4) we gave only an upper estimate for the hmit distribution of 
random walks analysing statistics of Markov chains on " locally free groups" . 

The consideration of problems dealing with the limit distributions of Markov chains 
on braid groups Bn requires examination of the "target space" of this group, i.e., the 
space where the random walk takes place. The structure of the target space is uniquely 
determined by the group relations and in the general case of the group -B„ > 4), is 
still unknown. 

In the present work we study the target space of the braid group Bn when n ^ 1, 
trying to develop a new "statistical approach" for words enumeration in this group. 

We should stress that our presentation offers a mathematical analysis which is far 
from rigorous, and ideas expressed here are mainly supported by numerical simulations. 
Moreover, we skip here some important but hard questions, like the problem of "words 
identity" in the braid group (deep advances concerning this subject can be found in 
recent work |S^). Our aim is to describe a constructive algorithm which, of course, has 
to be be justified and verified later. 

The structure of the paper is as follows: in the next section we give some necessary 
definitions concerning braid and "locally free" groups and describe the model under 
consideration; Section 3 is devoted to developing a "symbolic dynamics" method for 
words enumeration in the locally free group CJ-'n (for n 1); the target space of 
the braid group is studied in Section 4 by means of a statistical approach based on the 
concept of " locally free group with errors" . In this section we discuss also some additional 
links between this problems and the conventional random matrix theory, localization 
phenomena and statistics of systems with "quenched" disorder; we present here some 
speculations dealing with the possible relation of particular problems of random matrix 
theory to the theory of Dedekind function and modular groups. 

2 Basic Definitions and Statistical Model 

We first recall some points concerning the definitions of braid and "locally free" groups. 

Braid Group. The braid group Bn of "n-strings" has n—1 generators {(Xi, (72, • • • , cr„_i} 
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with the following commutation relations: 

(TjO-j+iCTi = CTi+iCTiO-j+i (1 < i < n - 1) 

aiCj = ajai (K - j| > 2) (1) 

aiCTj"^ = (Tj^Vj = e 

Let us mention that: 

- A word written in terms of "letters" — generators from the set {ai, . . . , ct„_i, erf ^, . . . , 
^n-i} gives a particular braid. Schematically the generators and may be repre- 
sented as follows: 




1 2 ••• i i + 1 ■ ■ ■ n - 1 n 




1 2 ••• z i + n 

- The length of the braid is the total number of letters used, while the minimal irreducible 
length hereafter referred to as "primitive word" is the shortest noncontractible length 
of a particular braid remaining after all possible group relations Eq.(|l]) are applied. 
Diagramatically, the braid can be represented as a set of crossed strings going from the 
top to the bottom after gluing the braid generators. 

- The closed braid appears after gluing the "upper" and the "lower" free ends of the 
braid on the cylinder. 

- Any braid corresponds to some knot or link. So, there is a possibility to use the braid 
group representation for the construction of topological invariants of knots and links, 
but the correspondence between braids and knots is not mutually single valued and each 
knot or link can be represented by an infinite series of different braids. 

Locally Free Group. The group CJ-'n{d) is called locally free if the generators, 
{/i, . . . , fn-i} obey the following commutation relations: 
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(a) Each pair {fj, fk) generates the free subgroup of the group CJ-'n{d) if \j — k\ < d; 

(b) fjfk = fkfj for \j -k\>d 

We will be concerned mostly with the case d = 2 for which we define CTn{2) = CTn- 

- The length of the word written in terms of letters {/i, . . . , /„_i, /f^, . . . , fn-\\ is the 
total number of generators used, and the "primitive word" is the shortest noncontractible 
length of a particular word after applying all relations of the group £JF„(2) (compare to 
the case of the braid group). The graphical representation of generators gi and gl^^ is 
also rather similar to that of braid group: 




1 2 ■■■ i i + n 




1 2 ■■■ % i + n 

It is easy to understand that the following geometrical identity is valid: 

^ ■ 'C^ 

% i+1 i i+1 

hence, it is unnecessary to distinguish between "left" and "right" operators /j. 

It can be seen that the only difference between the braid and locally free groups 
consists in the elimination of the Yang-Baxter relations (first line in Eq.(|l])). 

Statistical Model. Our aim is to calculate a specific "partition function", Ki(yU, ci), 
giving the number of all nonequivalent primitive words of length [i in the groups CJ^n+i {d) 
and Bn+i for n ^ 1. 
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Remark. To have a geometrical picture of the group CJ^n+i let us describe the recursion 
procedure of raising the graph (the "target space") associated with this group. 

Take the free group r„ with generators {/i, . . . , /„} where all /j {1 < i < n) do not 
commute. It is well known that the group r„ has the structure of a 2n-branching Cayley 
tree, C(r„) — see Fig.|l|a — where the number of distinct primitive words of length /i is 
equal to the number V^(/i) of vertices of the tree C(r„) lying at a distance of yU steps 
from the origin: 

Vn{li) = 2n{2n - 1)^"^ (2) 

The graph C{CJ-'n+i) corresponding to the group CJ-'n+i can be constructed from the 
graph C(r„) in accordance with the following recursion procedure: 

(a) Take the root vertex of the graph C(r„) and consider all vertices on the distance 
fi = 2 from it. Identify those vertices which correspond to the equivalent words in 
the group CJ-'n+i- (See example in Fig.|l|b). 

(b) Repeat this procedure taking all vertices at the distance /i = (1, 2, . . .) and "gluing" 
them at the distance fi + 2 according to the definition of the locally free group. 

By means of the procedure described, we raise a graph ("target space") corresponding 
to the locally free group CJ-'n+i- Now our main problem can be reformulated as follows: 
how many distinct vertices has the graph C [CJ-'n+i) at a distance of fi steps from the 
origin (for n ^ 1). 

In the next Section we give an exact answer to that question, which we use, in turn, 
as the basis for consideration of the much trickier case of the braid group. 

It is worthwhile to mention that the graph C{Sn+i) of the complete commutative 
group Sn+i (all generators of £n+i commute with each others) has the topology of the lat- 
tice embedded in R^'^ (see Fig.|l|c). Hence, the number of nonequivalent words V„'^°™™(/i) 
of length /i can be roughly estimated as the number of lattice points lying on the surface 
of the 2?T,-dimensional sphere, i.e. 

V~ifi) ~ const /i^" (3) 
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Comparing (j^) and (0) we get 



lim -lnK(yu) ^ ln(2n - 1) > 

n = const ^1 / , s 

< . (4) 

lim -lnV;"(^) = 

^ n = const 2> 1 

Naively we could expect that the behavior 

lim — InVnifi) = 

/i — > oo jJ, 
n = const 2> 1 

for "locally free" and braid groups remains unchanged (i.e. is the same as in the com- 
pletely commutative case) because for n ^ 1 we have of order ~ commutative 
relations and only of order ~ n noncommutative ones. However, it is proven for CJ-'n 
and conjectured for i?„ that 

lim — In = const > 

fl ^ oo H 
n = const ^ 1 

which clearly reflects the hyperbolic character of these groups. 



3 Exact Words Enumeration in Locally Free Groups 

We derive here an explicit expression for the number d) of all nonequivalent primi- 
tive words of length /i in the group CJ^n+i{d) (when d = 2 and n ^ 1) on the basis of the 
so-called "normal order" representation of words proposed by A.M. Vershik in and 



developed in [27, 28] which is reminiscent of the enumeration of "partially commutative 



monoids" known in combinatorics 



3.1 "Normal Order" Representation of Words 

Let us represent each primitive word Wp of length /i in the group CJ^n+i{d) in the 
normal order similar to the so-called "symbolic dynamics" appearing in the context 
of chaotic systems (see, for instance, 0) 

Wp = iuTHur ■ ■ ■ iur (5) 

where Y.i=i = ("^i T^OVi; l<s</i) and the sequence of generators /q,- in 
Eq.(^ for all distinct /q,. satisfies the following local rules [27|: 
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(i) If = /i, then fa,+, e {/2, /s, • • • /«}; 

(ii) lifa, = fkil<k<n - 1), then fa,+^ e {fk-d+i, • • • , /fc-i, /fc+i, • • • /n}; 
(in) If = /„, then G {/fc-d+i, . . . , 

These local rules could be represented diagramatically as follows: 



/i 




/l fs fn f2 fi fn fn-2 fn fn-l 

The rules (i)-(iii) give the prescription how to encode and enumerate all distinct 
primitive words in the group CTn+iid). If the sequence of generators in the primitive 
word Wp does not satisfy the rules (i)-(iii), we commute the generators in the word Wp 
until the normal order is restored. Hence, the normal order representation enables one 
to give the unique coding of all nonequivalent primitive words in the group CJ-'n+i{d). 

Example 1. Take an arbitrary primitive word of length // = 10 in the group CJ-'s+i{2): 
Wp 

— /s /s fs fi f2 fi fs fs fi fr 

= (/s)-^ (/s) (/s) ifi)-' (/2) if,) ifs? (/4) (/r) 

To represent the word Wp in the "normal order" we have to push all generators with 
smaller indices to the left when it is allowed by the commutation relations of the locally 
free group CJ-'q{2). We get: 

Wp = if,)-' ifsY {f,y if,)-' {f,f ifsf ifrY (7) 

(the "normal order" for this word is the sequence of used generators: {1, 3, 2, 5, 4, 8, 7}). 

To compute the number of different primitive words of length = 10 with the same 
normal order as in Eq.(|^), we have to sum up all the words like 

Wp = UiT' if^r' (/2)"^ (/s)™^ if^r^ Usr' u^r (§) 
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under the condition J2i=i = 10; 7^ V G [1, 7]. 



The calculation of the number of distinct primitive words, V^(/u), of the given length 
fi is now rather straightforward: 

Vn{fi,d) = J2Rn{s,d) ^ El^d-/^ (9) 

s=l {mi,...,ma} Lj=l 

where: 

• Rn{s,d) is the number of all distinct sequences of s generators taken from the set 
{/i! • • • ) /n} and satisfying the local rules (i)-(iii) 

• the second sum gives the number of all possible representations of the primitive path 
of length /i for the fixed sequence of generators — (see the example above); "prime" 
means that the sum runs over all mj 7^ for 1 < i < s; A is the Kronecker function: 
A = 1 for X = and A = for x ^ 1 

• special attention should be paid to the sequences built on the basis of one generator 
only, i.e. for primitive words of type Wp = {fkY ^ [l;"-] (see definition of Rn{s,d) 
below) . 

To get the partition function Rn{s,d) let us mention that the local rules (i)-(iii) 
define a generalized Markov chain with the states given by the nxn" incidence" matrix 
Mn{d), the rows and columns of which correspond to the generators /i, ...,/„ as it is 
shown below: 





fl 


/2 


h 


U 




fn-l 


fn 


/l 





1 


1 


1 




1 


1 




1 





1 


1 




1 


1 


fs 





1 





1 




1 


1 


h 








1 







1 


1 





































1 


fn 
















1 






(10) 



The matrix Mn{d) has a rather simple structure: above the diagonal we put everywhere 
"1" and below diagonal we have d—1 subdiagonals completely filled by "1"; in all other 
places we have "0" (the case with = 2 is shown in Eg . (|10|) ) . 
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The number of all distinct normally ordered sequences of words of length s with 
allowed commutation relations is given by the following partition function 



(11) 



where 



1 1 . . . 1 



and 



/ 1 M 

1 



> n 



V 1 / J 



(12) 



For s = 1 we have -R„(l, d) = VinVout = 2n as it should be. 

The remaining sum in Eq.(|) is independent on Rn{s,d), so its calculation is very 
simple (see also Appendix A for some generalizations): 



{mi,...,ms} 



.1=1 



2^ Cizl 



(13) 



Substituting Eq.(^ and Eq.(|TT|) into Eq.(g) we get 

vM = v^{fi,d) = X^2^c;:}/?„(s,rf) 



s=l 



(14) 



2vi„(2M„(d) + ir~W,^t = 2 Trace (2M„(rf) + J) 



where / is the identity matrix. 

Such a quantity is rather difficult to evaluate exactly. A reasonable approximation 



is to replace (|T^) by 



K(/i)=2E(2A. + l)^-^ 



(15) 



i=l 



where Aj are the eigenvalues of the matrix T„ which can be shown all to be real (see the 
later discussion) . In order to check the validity of approximation ( [I3| ) we have considered 
the case {n = 3, d = 2} where the exact value reads 

'15 -7v^' 



W^(^)(2-v/5r- 

whereas the approximation (|l5l) gives 



(2 + v^)^-i 



V3*(/i) = 2(2 - v^)^-^ + 2(2 + v^)^-^ + 2(-l)^-^ 

It can be seen that this approximation works reasonably well even for small values of fi. 
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The value VnifJ^, d) is growing exponentially fast with and the "speed" of this growth 
is clearly represented by the fraction 



q{d) 



Vn{lJ,,d) 



(16) 



which has the meaning of the effective coordinational number of the graph C{CJ^n 



In the next section we present calculations of the asymptotic expression of (|T5[) when 
n > 1. 

3.2 Calculation of Eigenvalues of Matrix M„(2) 



Consider the determinant 



,(A) = det (m„ - A J 



-A 1 1 
1 -A 1 
1 -A 



It satisfies the recursion relation 



a,(A) = -(A + l)a„_i(A) - (A + l)a„_2(A) 



with the boundary conditions 



ao(A) = 1 
ai(A) = —A 



For A > — 1 one may set 



which gives 



a„(A) = (A + l)^(-l)>„(A) 



ipn{\) = VA + lv9„_i(A) - V9„„2(A) 



(17) 



(18) 



(19) 



(20) 



(21) 



The general solution of (0) satisfying the previously defined boundary conditions 
19| ) is given in terms of Chebyshev's polynomials of second kind 



V9„(A) = W„+i(cosi9) 



(22) 
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where 
Therefore 

n — 1 

"^W = (-1)"(A + 1) — M„+i(cos<l) 

= (-ir(A+i)°i^^'°<"V>^ 

The last expression enables us to obtain all the eigenvalues of the matrix M„. In 
fact, it is convenient to distinguish them according to the parity in n: 

1. n = 2m + 1 

Ao = (— 1) (m such values); Afc=4cos^- fA;=[l,m+l]j (25) 

2. n = 2m 

Ao = (—1) {m such values); A^ = 4cos^ 1 (k = [1, ml) (26) 

2m + 2 ^ ^ 

Since in each case we have exactly n states, this exhaust the complete set of eigen- 
values, showing they all are real in the interval [—1,3]. One also recovers the result 
obtained earlier in |]27| , |2^ for the asymptotics of the highest eigenvalue of matrix M„ 
(in the limit 1): 



2 ^ 

Amax = 4 COS — - 1 

n + 2 



47r2 

3 - -t; (fc = 1) (27) 



n>l 



11? 



Now we are in position to compute the number of nonequivalent words V*{^) of the 
primitive length /i in the locally free group £JF„(2) for n ^ 1, n = const (see Eq.(|TB|)). 
Using the definition (|l^) and Eqs.(^5D-(^) we get for n = 2m + 1: 



= {n- l){-ir-' + 2 E I 8 cos^ - 1 (28) 



k=l 

Define (PnifJ') as follows: 



v:{fi) = {n-i){-ir-'+2^M 
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By means of Euler-Mac Laurin formula we may compute the asymptotic expression 

of ifnifl)- 



^ /"(Scos^x-lf-irfx-^T^-i 
2n Jo 2 

For for /i ^ 1, n=const^ 1 the last integral is evaluated by a saddle point approximation 
which yields 

Thus, for the number of nonequivalent words in the locally free group £JF„(2) we have 
the following limiting behavior: 

lim -lnV;*(/i) =ln7 (31) 

n = const ^ 1 

and Eq.(l^) gives 

q{d = 2) = 7 

Hence, the graph corresponding to the locally free group can be viewed as an effective 
tree with the branching number q = 7. 



(29) 



4 Approximate Statistical Approach for Words Enu- 
meration in Braid Group Bn 

The construction of an effective algorithm for enumeration of the words in the braid 
group Bn for n > 2 is of the most intriguing problems in group theory. 

In the present section we propose an approximate statistical approachioi enumeration 
of all distinct primitive words in the group i?„ for n ^ 1 which exploits some properties 
of locally free groups CJ^n considered above. 

The main idea is as follows. Let us deal with the sequences of words in the braid 
group Bn from the point of view of the locally free group CT^n "with errors". To be 
more specific let us start with the following example: 
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Example 2. Write a random word W in the group consisting of 8 letters. Let this 
word be for instance: 

W = ((7i)~^(74 ((75)~^ ((76)~^(75(Ti(76((72)~^ 

We reduce this word to the primitive one in two steps. 

1. On the first step we act in the same way as in the case of locally free group 

and push all generators with smaller indices to the left assuming that nearest neighbors 
do not commute at all. We get: 

W^reduced = {(^2)'^ CT^ {a^)'^ (f-e)"^ (75 ^^6 

^ V ' 

(76 ((75)-^ 

2. Now we can apply the Yang-Baxter relations to the triple {o'q)~^ (75 and obtain 
after the cancellation of ((75)"^ and (75 the primitive word 

Wp = (£72)^^4 CT6 (0-5)"^ 

The first step of the braid contracting procedure completely coincides with what wc 
did for the locally free group, while the second step we could regard (approximately, of 
course) as follows: 

Consider some pair, for instance, (erg) """^ (75. We commute it with the probability p. 
Such commutation we denote as an error. The probability to meet the generator Ui in 
the Markov chain with uniform distribution over the generators in the braid group is 
of order of p — Later on we consider more general case taking p as the variational 
parameter. 

So, let V^^^^'^{^) be the number of all primitive words of length /i in the braid group 
Bn. Our main idea is as follows: we would like to relate the quantity V^^'^^"^(/i) to the 
number of primitive words in the "group" CT^^^{2) averaged over the uniform distribu- 
tion of "errors" in commutation relations. 

It should be pointed out that the object CJ^^^{2) cannot be considered as a real 
group anymore. So everywhere below we understand under CJ^^^{2) just the ensemble 
of random words written in terms of the alphabet with specific commutation relations 
among the letters. 
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4.1 Statistics of Words with "Errors" in Locally Free Groups 



The methods of theoretical description of the systems with disorder are rather well 
developed, especially in regard to the investigation of spin glass models P3| . 

Central for these methods is the concept of self- averaging which can be explained as 
follows. Take some additive function F (the free energy, for instance) of some disordered 
system. The function F is the self-averaging quantity if the observed value. Fobs-, of 
any macroscopic sample of the system coincides with the value Fav averaged over the 
ensemble of disorder realizations: 

Fobs = {P) av 

The phenomenon of self-averaging takes place in the systems with sufficiently weak long- 
range correlations: in this case only F can be considered as a sum of contributions from 
different volume domains, containing statistically independent realizations of disorder 
(for more details see [ll8|). 

The central technical problem of systems with quenched disorder deals with the 
calculation of the free energy -F(/i) averaged over the randomly distributed quenched 
pattern. In our case we could associate the number of topologically different words with 
the partition function, hence the free energy would be Fi^jj) = — {\nV^^^{^)) and the 
"quenched pattern" is just the set of "errors" in commutation relations. 

A problem closely related to that mentioned above arises when averaging the correla- 
tion functions of some statistical system over the disorder. In this case the computations 
are based on finding the averaged density of states of some random matrix over the pre- 
scribed distribution of random entries. Below we show that the calculation of the mean 
value {V^"{fi)) belongs precisely to this class of problems. 

Conjecture 1 The number of nonequivalent primitive words, V^™*^(;u), of length fi in 
the braid group can be estimated in the limits n = const ^1, /i ^ 1 as follows: 

V^^^^\^^)^{V^{^^,d = 2,p)) (32) 

where Vn{fi,d = 2,p) is the number of all distinct primitive words of length fi with the 
"errors" in the commutation relations in the locally free group £JF„(2). We allow to 
commute the neighboring generators with the probability p and the averaging is performed 

16 



over the uniform probability distribution of "errors". The question concerning the choice 
of p is considered below. 



In support of our conjecture we bring the numerical computations performed in 
the work where we have constructed the (right-hand) random walk (the random 
word) on the group Qn = {'^^n^, Bn} with a uniform distribution over generators 
{gi, . . . , gn-i, gi'^ , ■ ■ ■ , gn-i} £ Gn- It mcaus that with the probability we have 
added the element ga^^ or g^^ to the given word of — 1 generators (letters) from the 
right-hand side. In ||28[ the following question has been raised: what is the averaged 
length of the primitive path (/i) for the A^-step random walk on the group 



In the Table 1 we show the results of numerical simulations carried on in p8| for the 
expectation value (fi) /N of the A^-step random walk on the "locally free structure with 
errors", CJ-''^^{2), and compare them to the same value for the A^-step random walk on 
the braid group Bn- 



Table 1. 





CT"^' withp= 1/5 


Bn 


n 






5 


0.55 


0.49 


10 


0.58 


0.56 


20 


0.59 


0.59 


50 


0.60 


0.61 


100 


0.60 


0.61 


200 


0.61 


0.61 



We have found asymptotically very good correspondence of the mean values (/i) /N 
for the braid group and the "locally free structure with the errors" for p = |- 

Of course, our conjecture is not exact and has a "mean-field" nature because we 
admit two subsequent generators {o'a^., o'ak+i) with nearest neighbor indices ak+i = afc±l 
to commute with the probability p regardless the value of the generator o"q,j,^2 in the 
sequence of letters in the word. Hence the problem of taking right value of p appears. 

Two limiting cases are trivial: 
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(i) For p = 1 we have a complete commutative group and the obvious inequahty is 
fulfilled 

C(/^)>(K(/i,ci = 2,p=l)) (33) 

(ii) For p = we return to locally free group CT21 for which we should have 

K''^"'^(/i)<K(/i,rf = 2,p = 0) (34) 



In the framework of the mean-field approximation and taking into account Eq.(32) 
we claim: 

(V;(/i,rf=2,p)) <V;^™<^(/x)<K(/i,rf=2,p = 0) (35) 
We show below that in the limit n 00, fi 00 the following equality takes place 



{Vr,{fi,d = 2,p)) 



Vn{fi,d = 2,p = 0) 



(36) 



The Main Statement. For the value V^^^^'^{fj,) we have the following asymptotic be- 
havior 



lim -InV^'^'^'^ifi) 

fj, ^ 00 jj, 

n ^ 00 



ln{Vn{fi,d = 2,p)) 



independent 
on p 



In 7 



(37) 



It is easy to understand that the number of nonequivalent primitive words VnifJ^, d = 
2,p) in the "locally free group with errors" can be calculated by means of averaging of 
Eq. (|n|) if we slightly change the matrix M„ replacing it by the random incidence matrix 
Mn: 



Mr{d = 2) 





fl 


/2 


/a 


/4 




fn-1 


fn 


fl 





1 


1 


1 




1 


1 


f2 


■^n—1 





1 


1 




1 


1 


/3 





Xn-2 





1 




1 


1 


/4 








•^n—3 







1 


1 


















fn-1 



















1 


fn 
















Xi 






(3^ 
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It is now parametrized by the random sequence of "0" or "1", i.e. 

= {x„_i,X„_2, • • ■,X2,Xi} (39) 

where 

Probfxo = 1) = 1 — » 

^ ' ' ^ (40) 

Prob(xj = 0) = p 

4.2 Density of States of Random Operator and Averaged Num- 
ber of Nonequivalent Words 

The determinant of the matrix M^"((i = 2) — A/ (see Eq.(|55D) satisfies the random 
recursion relation 

ak+i + ak{X + Xk) + ak-i{l + X)xk = (41) 
Introducing the Ricatti-like variable 

Pk = 

we arrive at the recursion relation 

Pk = -iX + Xk)-^-^^^ ke[0,n] (42) 

Pk-i 

with the boundary condition po = 

For the continuous sequence of {!}, i.e. {x*^"^} = {1111 ... 1} we have the 
following (nonrandom) transformation 

Pfc+i = -(A + l)(^l + -i^ (43) 

As soon as a zero appears in the random sequence, pk in ( ^2] ) is set to —A which coincides 
precisely with the initial value po- Since one returns back to the initial value, the process 
can be easily iterated for arbitrary random sequences {x^"-*} when ^ 1. Such a 
property of the map (^2]) is equivalent to the factorization of the determinant a„(A) of 
the random matrix M'^'^{d = 2) — A/: 

Example 2. Consider the sequence 

= ( 1 1 . .. 1 ... 1 ... 1) (44) 

h ™i h 
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The corresponding determinant a„(A) factorizes: 



an{X) = I[%W IlMr"' (45) 

{Ij} {ruj} 

It is worth pointing out that a recursion relation similar to (^31) appears also in the 
study of the binary product of random 2x2 matrices where one of the matrix is singular 
37| . Such a structure also occurs in the case of Ising model in a random magnetic 



field — see PSI. Following Derrida and Hilhorst [53], one can write down the invariant 



measure associated to (|42|) : 

oo oo 

v{p) = E p(i - p)'sip - Pk) + (1 - p) E p'^ip - Po) (46) 

k=l k=l 

The first (the second) term comes from the complete sequences of {1} ({0}) of arbitrary 
length k = 1,2,.... From the invariant measure one can compute Ina(A) which can 
be interpreted either as a Lyapunov exponent or as the free energy (depending on the 
physical context). So, we get: 

Ina(A) = lim lna„(A) = / V{p) In pdp 

L (47) 

= E^^'(i-#"'in«^(^) 

k=l 

Returning to our problem, we may use this expression to write down the averaged 
density of states 

1 /9 

P^^^ = -7^Imlna(A) = Ep'(1 -P)'"V.(A) (48) 
^ (^-^ k=i 

where Pn(A) = — Im lna„(A) is the density of states of a pure system (i.e. without 

TT 

randomness) of length n. From the density of states we can find the average number of 
words in the limit n oo 

Define {V*{p)) as follows: 

{V*{p)) = Jim = 2/(2A + ir-'^)dX (49) 

where the integration over A runs over the whole spectrum. Let us repeat once more 
that the function p(A) is the density of states of the random matrix Mn{d = 2) — A/ 
(ri — > oo averaged over the disordered pattern {x'-"^}. 
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The content of Eq.(|i8D is as follows. The density of states p(A) can be obtained by 
averaging the spectrum of the regular case "weighted" with associated sequences of {1}: 



Prob < 



X = (0 1 1 1 1 ■ ■ ■ 1 1 1 0) = p^jl - pY (50) 

complete set of {1} 



One should also add the contribution coming from zero's energy state corresponding 
to sequences of {0}. The resulting expression reads: 



p(A) = 5:^p2(l-p)2'"-i|m5(A + l) + ^^^<5(^A + l-4cos^^— ^ 

COS 



(51) 



m=0 

which may be rewritten as 



m=o L k=i \ 2m + 3^ 



rn + 21 



1 — Tl °° L 2 J / 7, 

p{\) = -^5{\ + 1) + Y.p\^-pT E 5 A + 1-4cos^^) (52) 

P n=0 k=l \ n + 6J 

where [x\ denotes the integer part of x. 



Using Eg . (1521) we may check that the function p(A) is properly normalized: 



p(A) = p(A) = 1 



-1 



Returning to (^) we get 



{V*{^^,p)) = 2 (]-^\ {-ly-' + 2f:p\l- pY S^^M (53) 



2-p 



n=0 



where 

Sn{f^) = j:Ucos'^-l\ (54) 



k=l 

(compare to (pSD ). 

In order to check the algebra we have computed (V^*(p)) for small values of p. One 
gets: 

{V*{1)) = {V*{2)) = 2 



which can be readily obtained through a direct calculation of lim 2Trace ( 2M„ + / 
We are however mainly interested in the limit fi —>■ oo. 
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Using (p9|)-([3T|) we may resum the series (|53| ) by isolating the contribution n < n 



and n > 11. After some algebra one obtains for < n + 1: 

Sn{^^) = (n + 1) E ci_, {-iY-^ci;\ - \ [i^-^ + (-1)^} 
p=l ^ 

which gives 



-1)'^ (55) 



T'^ hq < n < fi 
Snifi) ^^^^ ^ \ (56) 

where uq is some constant of order of unity. The Eq.(|55D is plotted in Fig.||a,b. 
Thus, we can rewrite Eqs.(|53|)-(plD as follows 

(r*(/i,p)) = 2 (^] i-ir-' +j2p\i-prsnM + f:/(i -p)"^„+2(/x) (s?) 



2-p 







The corresponding behavior of the function Q(p|yu) where 

Q(pM = i^^^^^^ (0 < p < 1) (58) 

is shown in Fig.^ for few fixed values fi = {10, 30, 150}. 

The plot in Fig.|^ enables us to come to the following conclusion. If the number 
of "errors" is small {p 0^), the volume of the group grows exponentially with the 
Lyapunov exponent In 7 (for /i oo). For the arbitrary number of "errors", p, the 
corresponding Lyapunov exponent approaches the same value In 7 for all p < 1 in the 
limit ^ ^ oo and exhibits a singular behavior just at the point p = 1 (which corresponds 
to the completely commutative group). 

The asymptotic expression ( p7D allows us to conclude that the limit behavior of the 
function V*{fi) is independent on p, €]0, 1[ and is the same as for the locally free 
group CJ-'n without any errors. This fact supports our conjecture (|37|) . 



It should be emphasized that these results are expected to hold only in the thermo- 
dynamic limit n —>■ oo. It would be more desirable to consider the limit in which the 
number of generators n is kept fixed and the length of the word, /i, is much larger than 
n. 
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4.3 Functional Equation, Continued Fractions and Invariant 
Measure 



The behavior of the spectral density of our model is very similar to the one encountered in 
the study of harmonic chains with binary random distribution of masses. This problem, 



which goes back to Dyson has been investigated by Domb et al and then thoroughly 
discussed by Nieuwenhuizen and Luck [Q. One considers a chain of oscillators where 
the masses can take two values: 

m with the probability \ — p 
M > m with the probability p 

In the limit M ^ oo the system breaks into islands, each of which consisting of n 
light masses surrounded by two infinite heavy masses. The probability of occurence of 
such an island is — p)". There is clearly a mapping to our model if one replaces the 
sequences of heavy and light masses by the sequences of "0" and " 1" . Many results may 
therefore be borrowed from the works 0, |3^ . In particular, by adapting the calculations 



of Nieuwenhuizen and Luck to our case one may rewrite the integrated density of states 



A 



^f(\) = I p{\')d\' 

in the form 



Ar(A) = l--^E(l-^>r^^) (59) 



where the relation between A and ^ is given in Eq. (|23|) . 

1 — p 



For A ^ — 1 one gets J\f(\) which corresponds to the contribution of the 

2 — p 

states A = — 1 at the bottom of the spectrum. 



At the upper edge of the spectrum (namely for A ^ 3^) one gets A/'(A) —>■ 1 which 
means that all the states are counted. Equation (|59|) shows that the behavior around 
A = 3 (corresponding to -t? = 0) is in fact dominated by the first term (n = 1) of the 
series. One has: 

P 



AfiX) ^ i--J—-{i-p) 



l-pV 

(60) 



- 1 P 



ln(l — p) 



V3^ 
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The behavior (|60| ) signals the appearence of Lifshits' singularity in the density of 
states. A more precise analysis shows that this result is in fact modulated by a periodic 
function |p5 |. 



Equation (|59D displays many interesting features. In particular, the function A/'(A) 
occurs in the mathematical literature as a generating function of the continued fraction 
expansion of — . 

V 

Let us briefly sketch this connection. Consider the continued fraction expansion 

I- ^ (-) 



Co + 



Cl + 



C2 + 



where all c„ are natural integers. Truncating this expansion at level n we get a rational 
number — which converges to — when n oo. 



A theorem of Bohmer |^6[ states that the generating function of the integer part of 

TT 

n=l 

is given by the continued fraction expansion 

z 1 



Giz) 



A2 + 

where 



AJz) 



/ 1\ 1n-2 

z) \Z, 



- 1 



and g„ is the denominator of the fraction — approximating the value — . 

In order to make connection with our problem it is sufficient to set 2; = 1 — p and 



express A/'(A) in terms of G{z). 

Equation (p9D shows that the invariant measure is a very singular object. However 
it satisfies a simple functional equation reminiscent of that which arises in the theory of 
automorphic forms. The equation can be derived either by a Dyson-Schmidt approach 
(see, for instance |T^) or just by looking at the explicit expression of 'P(p). 
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It is in fact simpler to work with the rescaled variable 

1 



which satisfies the recursion relation equivalent to Eq.(|4 



Zn = — (62) 

^n—l 



where n = vA+T. The transformation ( |62D may be obtained from the two matrices 
belonging to the group SL{2,R) 



( 



T 



1 /. _ s=\^ ~^ 
Oil \ 1 



The 517(2, R)-transformation TS = acts on z by the fractional linear transfor- 
mation 

azn^i + b 1 

= 7; — = ^ " 7 — • (^3) 

CZn-l -r Ct Zn-l 

The invariant measure of (0), which may be rewritten as 

00 

V{z) = E P(l - pT^ - [STTzo) (64) 

n=0 

can easily be shown to satisfy the fundamental equation 

V{z) = {cz + df V{Sfz) (65) 

up to some singular terms. 

By suitable rescahng it is in fact possible to absorb the prefactor 1/(1— p) and rewrite 
Eq.(|5D as 

V{gz) = {cz + d)-^V{z) (66) 

where g G SL(2,R). 

An analytical continuation of this expression into the complex z-plane would eventu- 
ally permit one to interpret V as an automorphic form. /^From the theory of automorphic 
functions it is well known that Eq.(^) is satisfied by the Dedekind modular function 

ri{z): 

^(^) = (^JT^^(^) M-^- = i) 
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Such objects although perfectly smooth in the upper half-plane Imz > display 
highly non-trivial fractal behavior on the boundary Imz = (see, for instance 0, [39(1). 

Another interesting connection which would be worth investigating is the fact that 
S and T generate the so-called Hecke group T{h) for h = 2cos^ (f? > 3 is integer). 
Surprisingly, these values of h coincide with a subset of the spectrum of the matrix M„ 
(see Eq.(ll)). 



5 Final Remarks 

5.1 The Geometrical View on the Word Enumeration Problem 

The number of primitive words in the locally free or braid groups allows a rather 
straightforward geometrical description. Namely, the matrix Mn{d = 2) can be re- 
garded as the transfer matrix for the model of a "biased Levy-flight" -like ("BLF"-like) 
one-dimensional random walk on the finite support. Actually, let us compute the sta- 
tistical sum of the process described below. Take n integers on the line: 1,2, ... ,n and 
consider the random walk when the walker can jump with equal probabilities from the 
vertex with the coordinate mi (1 < mi < n) to: 

(i) each vertex with the coordinate m2 (m2 G [mi + 1, n]); 

(ii) the vertex with the coordinate m2 = mi — 1 (i.e. one step back). 
The corresponding process is represented schematically in Fig.^. 

Analogously, we can associate the random operator M^^ with the transfer matrix of 
the generalized BLF-like random process which is described via the same rules (i) and 
(ii) but with additional requirement that the jump (ii) is blocked with probability p and 
allowed with probability 1 — p, independently of the position of the vertex. 

5.2 Conclusion 

We have proposed a statistical method for enumerating the primitive words in the braid 
group Bn based on the consideration of locally free groups with errors in commuta- 
tion relations. We brought arguments in support of the conjecture that the number of 
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long primitive words in the braid group is not sensitive to the precise local commuta- 
tion relations. We discussed the connection of the abovementioned problems with the 
conventional random operator theory, localization phenomena and statistics of systems 
with quenched disorder, and showed the connection between some particular problems 
of random matrix theory and the theory of automorphic functions. 

We believe that the problem of discovering the integrable models associated with the 
proposed locally free groups and developing the corresponding conformal field theory 
could help establish a bridge between the statistics of random walks on the noncom- 
mutative groups, spectral theory on multiconnected Riemann surfaces, and topological 
field theory. 

Acknowledgments 

We are very grateful to J. Desbois for elucidating for us many questions concerning the 
limiting behavior of random walks on locally free and braid groups; we would like to 
thank as well M. Tsypin for useful remarks and for help in the numerical confirmation 
of some of our conjectures. S.N. acknowledges the fruitful discussions with S. Fomin, L. 
Pastur, Ya. Sinai and A. Vershik on many aspects of the work. We highly appreciate 
the assistance of O.Martin in the final preparation of the paper and critical reading of 
the manuscript. 



27 



Appendix A 



The calculation of the number of distinct primitive words V^{fi) can be easily extended 
to the case of the group CJ-"^^^^ which is defined by the old relations: 

(a) Each pair {gj,gk) generates the free subgroup of the group ^JFP''™ if \j — k\ < 2; 

(b) gjgk = gkgj for \j -k\ > 2 

completed with the additionnal requirement^ {gtY = 1 which may be rewritten also as: 



gj = 1 |mod p\ Vj G [1, n] 



(Al) 



The calculation of the number of distinct primitive words, V^(/i), can be carried out 
using Eq.(P): 



J2 ~ ^^ 



(A2) 



vr'^i^^) = Y.Rnis,d) E ^ 

s=l {r?ii,...,r?is} Li=l 

where "double prime" means that the sum runs over all 7^ consistent with Eq.(A.l). 

Taking into account that the summation over all {mi, . . . , m^} obeying the condition 
(A.l) is independent on computation of Rn{s), we have to replace the Eq.(|13D by the 
following one 



TV 



{mi, . . . , rus} 
|mod p\ 



li=l 



2" c I ^"^ 

2TTi J \ 
(C) \™=i 



— { dz z-^-^+' ( 

(C) ^ 



where the contour C encloses the origin of the complex plane z 
We used the integral representation of Kronecker A-function 



(A3) 



A fx) 



dz 

2TTi J z^^^ 



x^O 

1 x = 



(AA) 



Let us consider three special cases: 
^Compare to the definition of the group of permutations. 
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1. p ^ oo. From Eq.(A.3) we get 

which coincides with Eq.(|l^). 

2. p = 2. From Eq.(A.3) we get 

3. p = 3. From Eq.(A.3) we get 

Substituting Eqs.(A.5)-(A.6) and Eq.(0) into Eq.® we may easily compute the 
corresponding values 
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Appendix B 



The computation of the averaged density of states p(A) of the random matrix M^" — AJ 
belongs to a class of problems with the "quenched disorder" (see the discussion in the 
section 4.1). Another class of problems closely related to the mentioned one arises when 
averaging the partition function (but not the free energy) over the disorder. Problems 
corresponding to the case of annealed disorder usually seem to be simpler from compu- 
tational point of view than that of quenched disorder but usually the thermodynamic 
behavior of systems with annealed disorder is less rich and less interesting than that of 
systems of quenched disorder. However in our case it would be very desirable to compute 



the value a„(A) where an{X) is determined by Eq.(^5|) in order to compare the values 
lna„(A) and lna„(A). 

Averaging a„(A) over the random distributions of {1} and {0} in the sequence 
(see (PD) and taking into account the factorization of corresponding determinant (^Sj), 
we get 

no s 



a„(A) 



-A) 



m,- — 1 



E E 

s=l {li...ls} {mi...ms} j=l 

where Y^' means that summation runs over the sequences {li...ls} and {mi 
obeying the conditions: 

s 

= n{l-p) 

3=1 



(5.1) 
(5.2) 



E 



np 



Recall that p is the fraction of zeros in the random matrix M^" and according to Eq.(B.2) 
the quantity {np) is always integer. 

Using the integral representation of the Kronecker A-function (see Eq.(A.3)), we 
may rewrite Eq.(B.2) in the form of a grand canonical distribution: 



an(A) = E(-Ar"^E n%+iw^ 

S = l Zi.../sj = l 



E h - ^(1 - P) 



' 27ri J 



E 

s=l L 



-A)-i^aKA); 



1=0 



(5.3) 



2TTi J 
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where the function b{X, z) plays a role of generating function for a; (A): 

oo 



Recall that 



ai{\) 



2{-l\ 



^(3-A)(l + A) 
1 



1 + A)'/2 gin + 2)6] 



sin 6* = --\/3 — A 



(S.4) 



cos 6* = -vT+A 

Collecting all equations together, we arrive at the following expression 
WAy= f-Ar^— ddzz^'^^-P^ A(l + z-^) + z-^ 



(5.5) 



Now we can use the definition of the Chebyshev's polynomials generating function: 

1 



l-2tx + ^2 



k=0 



where 



Unix) 



sin [{n + 1) arccosx] 



sin[arccos x\ 

After rather simple algebra we get the desired equation for the averaged determinant: 



a„(A) = (-1)"A"P- 



where 



t(i-p) 



Un{l-p){y) 



A3/2 



' A + 1 
A(A + i: 



1/2 



Un[l-p)-l{y) 



y 



2VA2^ 

As we can see, this result differs dramatically from the "quenched case" (see Eq, 



(5.6) 



(5.6a) 
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Figure Captions 



Fig.l. Graphs, corresponding to: (a) free group F^; (b) locally free group LTn+\\ (c) 
complete commutative group. In case of locally free group the vertices A and B should 
be glued because they represent one and the same word in group CJ^n+i- 

Fig. 2. Plot of the function In S'„(/i) in two regimes (Eq.(p3|)) 

Fig. 3. Plot of the function Q{p\n) for three fixed values of /x = {10, 30, 150} — see 
Eq.®. 

Fig. 4. Schematic representation of the process associated with the "biased Levy-flight" 
(BLF)-like random walk. 
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